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ABSTRACT 

We describe IRAC 3.6-8 fxm observations and ground-based near-IR 3HK S photometry from Mimir 
and 2MASS of the massive double cluster h & \ Persei complete to J=15.5 (M ~ 1.3M ). Within 25' 
of the cluster centers we detect ~ 11,000 sources with J< 15.5, ~ 7000 sources with [4.5] < 15, and 
~ 5000 sources with [8]< 14.5. In both clusters, the surface density profiles derived from the 2MASS 
data decline with distance from the cluster centers as expected for a bound cluster. Within 15' of the 
cluster centers, ~ 50% of the stars lie on a reddened ~ 13 Myr isochrone; at 15'-25' from the cluster 
centers, ~ 40% lie on this isochrone. Thus, the optical/2MASS color-magnitude diagrams indicate 
that h & x P er are accompanied by a halo population with roughly the same age and distance as 
the two dense clusters. The double cluster lacks any clear IR excess sources for J< 13.5 (~ 2.7M Q ). 
Therefore, disks around high-mass stars disperse prior to ~ 10 7 yr. At least 2 — 3% of the fainter cluster 
stars have strong IR excess at both [5.8] and [8]. About 4 — 8% of sources slightly more massive than 
the Sun (~ 1.4M Q ) have IR excesses at [8]. Combined with the lack of detectable excesses for brighter 
stars, this result suggests that disks around lower-mass stars have longer lifetimes. The IR excess 
population also appears to be larger at longer IRAC bands ([5.8], [8]) than at shorter IRAC/2MASS 
bands (K s , [4.5]), a result consistent with an inside-out clearing of disks. 

Subject headings: Galaxy: Open Clusters and Associations: Individual: NGC Number: NGC 869, 
Galaxy: Open Clusters and Associations: Individual: NGC Number: NGC 884, 
Stars: Circumstellar Matter, Infrared: Stars 



1. INTRODUCTION 

The evolution of circumstellar disks sets strong con- 
straints on the initial conditions for terrestrial and gas gi- 
ant planet formation. In young stars, the presence of cir- 
cumstellar disks is inferred by near-to-mid infrared (IR) 
emission by dust in excess of the output of stellar pho- 
tospheres. While the vast majority of stars in ~ 1 Myr 
clusters have disks, the frequency of disks declines on ~ 
1-10 Myr timescales (e.g. Hillenbrand 2005; Young et al. 
2004; Mamajek et al. 2004; Haisch et al. 2001). Be- 
cause circumstellar dust is the evolutionary precursor to 
> 1 km-sized planetesimals, the timescale for the disap- 
pearance of dust emission, the 'disk evolution timescale', 
is an important constraint for planet formation models 
(e.g. Strom et al. 1989; Hillenbrand 2005; Alibert et al. 
2005; Currie 2005). 

Recent observations of young (< 10 Myr) stars suggest 
that the evolution timescale may depend on stellar prop- 
erties and locations within the disk. K band (~ 2fj,m) 
excess from the inner (r < Q.lAU) disk is relatively un- 
common around both massive and very low mass stars 
(Hillenbrand et al. 1998). The - 5 Myr Upper Sco OB 
association also exhibits a mass-dependent frequency of 
disks: sources with IR excess at 8 /im and at 16 fj,m ap- 
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pear to be more prevalent around K/M and B/A stars 
than for F and G stars (Carpenter et al. 2006). Spa- 
tially resolved, mid-IR observations of ~ 3 — 5 Myr T 
Tauri stars (McCabe et al. 2006) show that many disks 
with strong mid-IR emission lack near-IR emission, pro- 
viding evidence for a location-dependent disk evolution. 
The drop in near-IR circumstellar dust emission may 
be explained by grain growth, which is probably faster 
in the inner disk regions due to higher midplane densi- 
ties and orbital frequencies (e.g. Dullemond & Dominik 
2005). Circumstellar gas dispersal by UV photoevapora- 
tion (Clarke et al. 2001; Alexander et al. 2006) and plan- 
etesimal formation by gravitational instability (Youdin 
& Shu 2002; Youdin & Chiang 2004; Currie 2005) may 
also result in a location-dependent evolution and, thus, 
a wavelength-dependent timescale for the disappearance 
of near-IR emission. 

Multiwavelength near-to-mid IR observations of 
evolved (> 10 Myr), massive (> 1000 members) clusters 
are required to measure the disk evolution timescale as 
a function of stellar mass and disk properties. The dou- 
ble cluster h & % Persei, the most massive, evolved open 
cluster within 3 kpc, provides a rich laboratory to study 
disk evolution in more detail. Starting with the initial 
study of Oosterhoff (1937), there has been much debate 
concerning the age, distance, and stellar content of the 
clusters (e.g. Borgman & Blaauw 1964; Wildey 1964; 
Schild 1967; Crawford et al. 1970; Vogt 1971; Tapia et 
al. 1984; Marco & Benabeau 2001; Capilla & Fabre- 
gat 2002). Recent work (Keller et al. 2001; Slesnick et 
al. 2002, hereafter S02; Bragg & Kenyon 2005, hereafter 
BK05) has converged on a nearly identical age for both 
clusters, ~ 12-13 Myr, a common distance of ~ 2.34 kpc 
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(a distance modulus ~ 11.85), a low extinction of E(B- 
V)^0.5 uniform across the clusters (BK05), and an ini- 
tial mass function (IMF) of massive stars consistent with 
results for other young, massive clusters (Massey 2003; 
Bragg 2004). Although the clusters have similar masses 
of > 4000 M Q (h Per) and > 3000 M Q ( X Per), the in- 
ternal dynamical structure of the clusters may be very 
different (BK05), with some evidence for mass segrega- 
tion in both clusters (S02, BK05). 

While numerous optical photometric and spectroscopic 
studies of h & \ Persei exist, the cluster has yet to be 
explored in detail in the near-to-mid IR. Previous IR ob- 
servations have been shallow, J~K < 11-13 (Tapia et 
al. 1984; Bragg & Kenyon 2002, respectively), and con- 
centrated on understanding the population of Be stars. 
While some constraints have been placed on the stellar 
population of h & \ P er i the population of circumstellar 
disks remains unprobed. Thus, its potential to inform 
disk evolution theories remains untapped. 

In this paper, we describe the analysis of the first 
deep IR survey of h & \ Persei. The survey combines 
2MASS All-Sky Survey (2MASS) and Infrared Array 
Camera (IRAC) data covering the entire double cluster 
with deeper near-IR data of selected cluster regions from 
the Mimir camera at Lowell Observatory. With approxi- 
mate completeness limits in JHK S about 3-5 magnitudes 
deeper than BK02 and Tapia et al. (1984), respectively, 
the 2MASS survey reaches stars with masses smaller than 
the limits of all optical surveys except for Keller et al. 
(2001) with less interference from the large line-of-sight 
extinction. When combined with longer wavelength (3.6, 
4.5, 5.8, & 8 H Spitzer/IRAC data, 2MASS and deeper 
near-IR surveys allow probes of disk evolution for stars 
with ages ~ 10-15 Myr and masses > 1.3 Mq. From our 
analysis we hope to provide valuable input for circum- 
stellar disk and massive cluster evolution models. 

Our results provide good evidence for two populations 
of ~ 13 Myr old stars in the direction of h & \ Per. In 
both clusters, the sky surface density declines approxi- 
mately inversely with distance at < 20' from the cluster 
centers, which is consistent with results derived from op- 
tical data (Bragg & Kenyon 2005). At distances beyond 
~ 20' from the cluster centers, the surface density is con- 
sistent with the background density. CMDs indicate that 
nearly half of the stars within 15' of the cluster centers 
lie on a reddened 13 Myr isochrone. In regions of lower 
surface density 15 '-25' from the cluster centers, however, 
roughly 40% of the stars lie on the same isochrone. Thus, 
the double cluster may be accompanied by an extensive 
halo population, much like the Orion star-forming region, 
that has roughly the same age and distance as hone fide 
cluster members. 

While the majority of sources have photosphcric IRAC 
colors through 8/xm, the cluster and halo population as- 
sociated with h & x Persei have a small IR excess pop- 
ulation. Disks around stars with J < 13.5 (~ 2.7 Mq) 
are extremely rare: disks around massive stars disperse 
by ~ 10 7 years. The IR excess population is larger for 
fainter sources down to J=15.5; disks around lower mass 
stars may have longer lifetimes. The IR excess popu- 
lation also grows progressively larger with longer IRAC 
wavelengths, a result consistent with an inside-out clear- 
ing of protoplanetary disks. 

We begin with a description of the data and reduc- 



tion in Section 2 and discuss the analysis of the near-IR 
data in Sections 3 and 4. In Section 5, we develop our 
technique to identify IR excess souces and use this ap- 
proach to make several estimates of the fraction of IR 
excess sources in the cluster. We conclude with a brief 
summary in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Near-Infrared Ground-Based J,H, & K s data from 

2MASS 

2MASS is an all-sky survey with uniform, complete 
photometry at J, H, and K s (Skrutskie et al. 2006). The 
survey has 10cr sensitivity for point sources with J > 15.8, 
H > 15.1, and K s > 14.3. For sources with 10<r detec- 
tions, the survey is more than 99% complete and more 
than 99.95% reliable, with la astrometric accuracies of at 
least 0.1 arcsec relative to the Hipparcos reference frame 
for sources with K s < 14. Although the 2MASS survey 
used cameras with 2 arcsec pixels, multiple observations 
of each sky position yields images with a nominal reso- 
lution of 1 arcsec. 

We downloaded data for h & \ Persei from the 2MASS 
point source catalog (PSC) using the IRSA web inter- 
face at IPAC 6 . The catalog includes ~ 31,000 sources 
with 5a detections within 1 deg of aj2ooo = 2 ,l 20 m 2.9 s 
<Sj2ooo = 57°5'41.7" (1=135°, b=-3.7°), which is a point 
midway between the two clusters. A skymap from the 
2MASS PSC clearly shows the two clusters as dense 
concentrations of stars (Figure 1). There are ~ 11,000 
sources within 25' of either cluster center. 

Figure 2 shows the magnitude and error distributions 
for J and K s using stars within 25' of the cluster centers. 
At J (K s ), the number counts monotonically increase 
to J = 15.5 (K s = 14.8) and then turn over. Both of 
these peaks are somewhat brighter than the nominal 10cr 
sensitivity limits, suggesting that confusion of sources in 
the galactic plane or the cluster centers might limit the 
source counts. However, both clusters show secondary 
peaks at J ss 14.5 and at K s 14. We show below that 
stars within the two clusters produce these peaks. 

2.2. Near-Infrared Ground-Based J, H, & K s data from 

Mimir 

The relatively low spatial resolution and the shallow 
peaks in the JHK S magnitude distributions of the 2MASS 
data prompted us to acquire additional near-IR photom- 
etry^ We used Mi mir, a m ultifunction IR instrument 
( |http: / /people.bu.edu/ clemens/mimir/[ ) at the f/17 fo- 
cus of the 1.8m Perkins telescope at Lowell Observatory 
(Table 1). Mimir uses a Mauna Kea JHK filter set and 
covers a 3'x3' field with 0.18" pixels in this configura- 
tion. The total on-cluster coverage was ~ 144 square arc- 
minutes (Figure 1). A series of flats, darks, and biases 
was taken at the beginning and end of each night. The 
telescope was dithered by 30" for each pointing where a 
series of three 10-second exposures were taken. In most 
cases the pointings were too crowded to construct sky 
frames from median filtering. Therefore, we took a se- 
ries of three 10-second dithered pairs of exposures for sky 
frames 7' in right ascension and declination away from 
each dithered pair of on-cluster pointings. 

6 http://irsa.ipac.caltech.edu/ 
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We followed a standard image processing procedure. 
First, we subtracted a dark frame from each image and 
divided by a dark-subtracted and normalized flat field. 
Next, we subtracted each object frame by the appropri- 
ate median-filtered sky frame. We used a custom IDL 
procedure to interpolate over bad pixels as well as im- 
age/detector artifacts unique to each frame. We used 
another custom IDL procedure (similar to 'imcombinc' 
in IRAF) to stack 10-second frames together. Up to 18 
individual frames were stacked together for each pointing 
yielding a total integration time of ~ 3 minutes each in a 
given field. We obtained longer integrations of ~ 7 min- 
utes for two of the fields, one each adjacent to the h & \ 
Per centers, typically improving the completeness limits 
by ~ 0.3 magnitudes. Each stacked image was visually 
inspected for errors in offset computation and residual 
image artifacts. 

We identified sources and extracted aperture photome- 
try using SourceExtractor (SE; Bertin & Arnouts 1996). 
Each potential source was convolved with a 9x9 gaus- 
sian filter (5-pixel FWHM) and sources on image edges 
were removed to eliminate spurious detections. We typ- 
ically used a 12 pixel diameter for aperture photome- 
try in all bands with a tendency to use smaller aper- 
tures for K data (due to smaller point-spread functions) 
and subtracted the background from a filtered, global 
background map. Sources from each dithered pair were 
matched to generate a final list of detections in each 
band. 

To derive an absolute calibration, we matched Mimir 
sources with 30-40 2MASS zero-color sources on each 
frame. Once we established the zero-point, these sources 
had < 1% offset as a function of JHK S . For redder 
sources, we measured a reliable offset between K(Mimir) 
and K S (2MASS) of ~ 3% on sources with J-K - 1. Ap- 
plying the conversions from Carpenter et al. (2001) to 
correct Mimir colors to 2MASS colors eliminated the off- 
set. For the ensemble of Mimir sources, we estimate ab- 
solute uncertainties of < 0.025 mag relative to 2MASS 
in all bands. 

To derive coordinates for the Mimir sources, we re- 
lied on accurate 2MASS astrometry. Although tests with 
standard packages, such as WCStools, showed that the 
Mimir fields were well-oriented north-south and showed 
little distortion and rotation, the small Mimir fields made 
robust astrometric solutions difficult. Because our main 
goals for astrometry were matching 2MASS sources and 
estimating incompleteness, we computed J2000 coordi- 
nates for Mimir sources by matching several 2MASS 
sources per field and deriving relative coordinates using 
the known pixel scale. Comparisons of all Mimir sources 
with 2MASS counterparts yields an average positional 
offset of 0.8" ± 0.4", with little evidence for systematic 
offsets as a function of position on the Mimir detector. 
Although better coordinate accuracy might be possible 
with complete astrometric solutions for all Mimir fields, 
this positional uncertainty is sufficient for robust source 
matching even in the centers of each cluster. 

Figure 3 shows the J and K s magnitude distributions 
for the Mimir data. These data reach ~ 0.5 mag deeper 
than the 2MASS data, with clear peaks in the counts at 
J = 16 and at K s = 15.5. Of the ~ 1000 Mimir detec- 
tions, ~ 650 sources with J < 16 have 5<r detections in 
all three bands. Multiple observations of ~ 50 sources 



verified these uncertainties. At the 2MASS lOc sensitiv- 
ity limit, we recover all 2MASS sources. At this limit, 
the typical magnitude difference was <5J ~ 0.1 mag and 
SK S w 0.15 mag, which is consistent with the expected 
error distribution. 

Figure 3 shows the same secondary peaks as the 
2MASS data. For a cluster age of 13 Myr, a distance 
modulus of 11.85, the Siess et al. (2000) isochrones, 
and the color conversion table from Kenyon & Hartmann 
(1995), the peaks at J = 14.5 and K s = 14.25 correspond 
to stars with masses M — 1.6M©. With these assump- 
tions, stars at the J = 16 limit have masses of 1.1—1.2 
Mq. In Section 3, we restrict our analysis of the Mimir 
data to those 650 sources detected at all bands. 

2.3. IRAC 3.6-8 fim data 

We obtained observations of h & % Persei on Jan- 
uary 18, 2004 with IRAC (Fazio et al. 2004) on the 
Spitzer Space Telescope. The IRAC survey covers about 
~ 0.75 deg 2 centered on 02000 = 02:20:29.166, 820m = 
+57:12:27.54 (Figure 1). There is a ~ 7 arcmin offset 
between the channel 1/3 and channel 2/4 mosaics. We 
used the 12s high dynamic range mode to obtain two 
frames in each position, with 0.4s and 10.4s integration 
times. The observation of each field was repeated with 
a small offset, which allowed 20.8s integration time for 
each position. We identified image artifacts and cosmic 
ray hits by comparing the two observations. 

The data were taken during a period of above-normal 
solar activity. To remove cosmic rays we evaluated two 
approaches. First, we took the conventional steps of 
mosaicing and extracting sources, relying on the cos- 
mic ray circumvention features of the software. We used 
PhotVis (v. 1.09) for source finding and aperture pho- 
tometry (see Gutermuth et al. 2004). The radii of 
the source aperture, and the inner and outer sky an- 
nuli were 2.4, 2.4, and 7.2 arc-seconds respectively. In 
the second approach, we used a custom IDL routine de- 
veloped by T. C. to extract sources individually from 
the BCD data prior to mosaicing, applied array-location- 
dependent photometric corrections (see Quijada et al. 
2004; http:/ /ssc. spitzer. caltech.edu/irac/locationcolor), 
removed sources near image edges, and interpolated over 
pixels flagged as cosmic ray hits/image artifacts in the 
first method. This approach degraded the detection limit 
but allowed better control of the cosmic ray artifacts than 
by mosaicing alone. Source Extractor aperture photom- 
etry was performed using a 4 pixel diameter aperture; 
the background count level and rms was computed from 
the filtered, global background pixel map. 

We compared the results of the two methods and exam- 
ined the raw data to determine the best method for each 
band. This comparison included examining the number 
of sources with colors bluer than Rayleigh Jeans and the 
scatter in color-color plots, both taken as indicators of 
cosmic ray effects. The source extraction prior to mo- 
saicing provided the most reliable results in bands 1 and 
2 ([3.6] and [4.5]), while mosaicing first was better in 
bands 3 and 4 ([5.8] and [8]). The slightly larger PSFs 
in the latter two bands provided enough pixels for the 
circumvention software to distinguish cosmic rays from 
real point sources. Since the intrinsic ratio of signal to 
noise is higher at [3.6] and [4.5] than [5.8] and [8], the 
degradation from extracting sources prior to mosaicing 
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was not serious for the shorter wavelength bands. 

We calibrated the photometry using large aperture 
measurements of standard stars obtained during h & 
X Per observations, applying an aperture correction for 
each channel to account for the difference between the 
aperture sizes used for standard and h & \ P er photom- 
etry (see Reach et al. 2005 for calibration details). The 
brightest sources (< 10th magnitude) are saturated even 
in the short exposure frames. Figure 4 shows the mosaic 
image of the [3.6] channel at low contrast. 

We bandmerged the data using a sub-arcsecond match- 
ing radius to minimize the contamination of any residual 
cosmic ray hits mistakenly identified as 'sources' by the 
source extraction algorithms of SE and Phot Vis. 

To remove sources that are likely AGN or galaxies with 
aromatic emission, we rely on published optical data, 
our 3KK S data, and IRAC colors. Requiring 2MASS, 
Mimir, or optical counterparts eliminated most highly 
reddened extragalactic sources. In the IRAC bands, typ- 
ical field AGN have [4.5] > 14 and [4.5]- [8] > 1.25 (Guter- 
muth, unpublished); galaxies with aromatic emission lie 
to the right of a line from ([3.6]-[5.8],[4.5]-[8]) - (0,1) to 

- (1.5,3) with [4.5]-[8] > 1 as well as ([4.5]-[5.8],[5.8]-[8]) 

- (0,1) to - (1,2.25) with [5.8]-[8] > 1. We identify very 
few sources with optical/near-IR counterparts and [4.5]- 
[5.8] > 1 or [5.8]-[8] > 1. Thus, extragalactic sources 
have negligible impact on our analysis. 

Figure 5 shows magnitude distributions for the four 
IRAC bands and error distributions for [4.5] and [8]. At 
[3.6] and [4.5], the distributions have a monotonic rise 
to [3.6] w [4.5] - 14, a broad plateau at [3.6] w [4.5] 

- 14-15, and a sharp drop at [3.6] w [4.5] > 15. The 
longer wavelength IRAC bands have steeper magnitude 
distributions and brighter magnitudes with peak count 
levels ([5.8] - 14.25 and [8] - 13.75). The errors shown 
in the IRAC bands demonstrate that the average pho- 
tometric uncertainty is well below 0.2 (~ 5<r) at least 
through [4.5]=15 and [8]=14.5. 

As with the 2MASS and Mimir data, the IRAC data 
show evidence for a secondary peak in the magnitude 
distributions. This peak lies at [3.6] w [4.5] ~ 14 and [5.8] 
<~ 13.75. Adopting the Siess et al (2000) isochrone and 
the Kenyon & Hartmann (1995) colors for main sequence 
stars, the peaks at JHK S and [3.6]-[5.8] are consistent 
with cluster stars having masses ~ 1.3 M Q at a distance 
modulus w 11.85. 

Figure 6 shows our estimate of sample completeness, 
where we plot the ratio of sources detected in both IRAC 
and 2MASS to those detected in 2MASS for each IRAC 
band as a function of 2MASS J magnitude (for J< 15.0). 
The sample oscillates between 90 and 100% complete in 
each IRAC band at magnitudes brighter than expected 
based on where the source counts peak in each band. 
This discrepancy, most noticeable in the [4.5] band, is 
due to the very small matching radius employed in com- 
bining the 2MASS and IRAC data sets, not photometric 
errors. If the matching radius is increased, the complete- 
ness level in the IRAC bands increases at the expense 
of larger cosmic ray contamination. Because our goal 
is to estimate the fraction of sources with small IR ex- 
cesses, we prefer to analyze a smaller sample with more 
robust colors. When the if s -[IRAC] colors are analyzed 
in Section 5, we consider the effect that uneven sample 
completeness has on our estimates of the IR excess pop- 



ulation. Table 2 lists our photometry from 2MASS and 
IRAC. 

3. GROUND-BASED mK s DATA ANALYSIS 

To make reliable estimates for the fraction of cluster 
members with IR excesses, we must (i) derive the frac- 
tion of stars in either cluster, (ii) derive robust crite- 
ria for defining an IR excess, and (iii) combine the two 
criteria into clear estimates for the excess fraction. In 
this section, we consider two methods for estimating the 
fraction of stars within each cluster, number counts and 
model isochrones. The first method measures the num- 
ber of stars in the clusters relative to the stars in the 
background population. Model isochrones allow us to 
estimate the fraction of stars with magnitudes and col- 
ors that are consistent with the magnitudes and colors 
of stars on a 13 Myr isochrone (S02). Both approaches 
yield similar results. 

With probabilities for cluster membership established, 
we then consider whether any cluster stars have near-IR 
excesses. These data demonstrate that the vast majority 
of sources brighter than the completeness limits have col- 
ors consistent with photospheric colors at 1-2 /im. Thus, 
we find no evidence for near-IR excess sources in h & \ 
Per for J < 15.5, which correspond to cluster stars with 
masses > 1.3 Mq. 

3.1. 2MASS archival JHK S data: density distribution 
and cluster membership 

3.1.1. Density Distribution of Sources on the Sky 

To investigate spatial inhomogeneities and other struc- 
ture, we computed the projected sky surface density of 
stars, a standard approach for deriving the properties of 
star clusters (e.g. Binney & Tremaine 1987). We re- 
strict our analysis here and in Section 3.1.2 to sources 
brighter than J=15.5. To derive the surface density of 
stars (top panel of Figure 7), we counted stars in 1.5 x 
1.5 arcmin bins, displaying the density in 10% increments 
from 0-90% of the peak surface density of 15 arcmin -2 in 
the center of h Per. This surface density is higher than 
the ~ 7-8 arcmin -2 peak found by BK05, who were re- 
stricted to sources with spectral types earlier than A5. 
Based on J magnitude to spectral type conversions us- 
ing the Siess et al. (2000) isochrones and preliminary 
MMT/Hectospec spectroscopy (Currie et al. in prep), 
our population includes sources with spectral types ear- 
lier than <~ K0. 

Aside from the strong peaks, each cluster has consider- 
able small-scale structure. BK05 report a kidney-shaped 
isodensity contour in \ P er j a fairly symmetric inner core 
surrounded by a rectangular isodensity contour in h Per, 
and other structures ~ 5-10 arcmin from the cluster cen- 
ters. With the deeper 2MASS data, these structures are 
clearly visible at ~ 40 — 60% of the peak density and well 
above the median background level of ~ 2.7 arcmin -2 . 
The mean level is <~ 2.8 arcmin -2 and the background 
noise fluctuates by cr~0.65 arcmin -2 . 

Other asymmetric structures are also apparent ~ 10-15 
arcmin from the cluster centers, which may suggest that 
both clusters are asymmetric on large scales, ~ 10 pc at 
a distance of 2.34 kpc. However, features > 15'-25' from 
the cluster centers have amplitudes comparable to the 
noise level and thus might not be real. Deeper near-IR 



h & x Persei 



5 



data could verify the existence of lower-amplitude struc- 
ture. 

The bottom panel of Figure 7 shows the radial surface 
density plots derived from the 2D map, using l'-wide 
half-annuli (facing away from midway point of h & \ 
Per). In both clusters, the surface density drops rapidly 
from ~ 2.5 to 10 arcmin and then slowly merges into the 
apparent background level ~ 30 arcmin from the cluster 
centers. The surface density reaches the median back- 
ground density level of 2.7 arcmin -2 (er=0.278 arcmin -2 ; 
J < 15.5) <~ 20-25 arcmin away from the h& \ Per cen " 
ters at 2H8 m 56.4 s , 57°8'25" and 2' l 22 m 4.3 s , 57°8'35" 
respectively (BK05). Although the clusters clearly are 
not symmetric, we can derive a reasonable estimate of 
the fraction of stars within the clusters by integrating 
the surface density above this background level. This 
approach yields a cluster population of ~ 2000 stars, ~ 
35% of the total population within 15 arcmin of the clus- 
ter centers. 

Although the median background level is reached ~ 
15 — 20' there is a small gradient in the radial surface 
density profile from ^10 — 35 arcmin (e.g. to the limit 
of our sample), though no clear general gradient across 
the field along constant right ascension or declination. 
There are also is at least one region beyond 10' from ei- 
ther cluster center (from ~ 12 to 20' away from % Per) 
that may extend above its surroundings background level 
on scales larger than size scale of the background fluctua- 
tions. This behavior motivates us to investigate the stel- 
lar content within the 'background' (beyond <~ 15' from 
the cluster centers). Thus, we now compare the color- 
magnitude (CMD) diagrams of the cluster-dominated 
and low-density regions. 

3.1.2. Possible Evidence for a Halo Population o/<~ 10 7 
year old sources in the vicinity of h & \ Persei 

With no evidence for a significant age spread (S02) 
or patchy extinction (BK05) across the clusters, we can 
use the cluster's isochrone to constrain cluster member- 
ship. We construct V/V-J and J/J-H color-magnitude 
diagrams, using data from S02 and 2MASS. For both 
color-magnitude diagrams, we adopt an age of 13 Myr for 
the cluster (S02) using theoretical isochrones from Siess 
et al. (2000). For the J/J-H diagram we also used pre- 
main sequence tracks from Bernasconi et al. (1996) for 
comparison, which show excellent agreement. To com- 
pute the cluster reddening, S02 restricted their analysis 
to the cluster nuclei and found a median reddening of 
E(B-V) <~ 0.56. Using a larger spatial sample of sources, 
BK05 computed a slightly smaller E(B-V) of - 0.52. 
Since the larger spatial sample of BK05 is more similar 
to ours, we consider this estimate more appropriate for 
our sample and adopt it. This choice has negligible bear- 
ing on our results since the reddening is low and yields 
an age within la of S02's value (Bragg 2004). Convert- 
ing from optical to infrared extinction via Bessel & Brett 
(1988) yields E(J-H) - 0.185 for E(B-V) - 0.52. 

If the clusters have a small spread in age, the frac- 
tion of stars identified as being on the 13 Myr isochrone 
should be similar to the ~ 35% estimate derived from 
the number counts in Section 3.1.1. For objects with J 
< 14.5, the V/V-J diagram provides an efficient method 
to test for cluster membership. The S02 V data have 
corresponding 2MASS J limits of <~ 14.5; the diagonal 



slope (see Figure 8, top panels) of the isochrone on V/V- 
J allows good detection of blue, early spectral type (typ- 
ically) background, non-member sources and very red, 
foreground non members. For objects with J > 14.5, the 
V data become incomplete. However, as shown in Figure 
8 (bottom panels), the J/J-H CMD isochrone is sensitive 
to membership in this magnitude range and thus pro- 
vides a good substitute for the optical CMD. The two 
CMDs are used together to identify sources of the same 
age and distance of h & \ P er brighter than J = 15.5. We 
divide our sample into two main populations: an 'on clus- 
ter' population corresponding to sources within 15' of the 
cluster centers and a low-density population for sources 
between 15' and 25' away. We chose 15' for the first pop- 
ulation because it corresponds to ~ 3-4 cluster core radii 
(BK05). Based on our density distribution analysis, this 
distance corresponds to the point where the typical stel- 
lar density begins to approach the median background 
level. The 25' outer radius for low-density regions is cho- 
sen primarily because it fully samples our entire IRAC 
coverage. According to the modified Hubble law distri- 
bution (p ~ r -2 ) of cluster members, the cluster density 
<~ 4 core radii away should be ~ 6% of the peak density. 
If the number of stars in a circle of diameter 1 arcmin at 
the center is N, an annulus of 1 arcmin width at 15 ar- 
cmin should have roughly 4-5 N stars in the cluster and 
roughly 25-30 times that many in the 'background'. The 
fraction of stars in the cluster beyond 15' to 25' should be 
over twice as small (p ~ 0.01-0.03p pea fe). Therefore, if h 
& x Persei sit in a background of non-related field stars, 
the low-density regions beyond 15' will be dominated by 
sources of different ages and distances than sources in h 
& x Per. The density of sources tracking the isochrone 
should be far larger within 15' than outside 15'. 

The V/V-J diagram for sources within 15' of the h 
& x Persei cluster centers (Figure 8, top left panel) 
clearly shows a distribution of sources tracking the 13 
Myr isochrone for V=12-16. The isochrone slightly bends 
at V ~ 15-15.5. Foreground sources arc clearly visi- 
ble above the isochrone. Figure 8 (bottom left panel) 
shows the J/J-H CMD diagram for all sources within 
15'. Sources with J < 14.5 track the reddened isochrone 
well except for some sources with J-H > 0.5 that are 
probably foreground M dwarfs or background M super- 
giants. The isochrone 'bends' horizontally at J ~ 14.5. 
Almost all sources to J = 15-15.5 fall along the 13 Myr 
isochrone or are redder. Although many sources with J > 
15.5 may follow the isochrone, the larger errors (a > 0.1) 
and increasing population of blue sources make it hard 
to measure the cluster population. Thus, we restrict our 
analysis to stars with J < 15.5. 

The top-right panel of Figure 8 shows the V/V-J di- 
agram for sources both 15'-25' from each cluster center 
and greater than 15' from both centers. While sources 
appear quite scattered around the isochrone for V < 14, 
there is a significant population of sources between the 
two isochrones (for single and binary stars) for V = 14 
- 16.5. There are also many sources between the cor- 
responding isochrones in the J, J-H CMD (Figure 8, 
bottom right panel). From analysis of the projected 
radial density distribution, regions > 15' from either 
center should contain mostly background or foreground, 
non-member sources. The fact that the fainter of these 
sources loosely follow the 13Myr isochrone for the cluster 
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suggests that there may be a larger region of enhanced 
formation of low mass stars in the same general direction. 
This behavior would be similar to the situation with the 
Orion cluster, for example. 

We now quantify the fraction of sources in the h & 
X Persei 2MASS coverage that are consistent with clus- 
ter membership by counting the number of stars that are, 
within photometric errors, lying on the V/V-J and J/J-H 
isochrones and the number lying outside them. We re- 
quire that sources be within 0.3 magnitudes (~ 5a added 
in quadrature) of the isochrone to count as a member. 
This procedure neglects colors expected for unresolved 
binaries, which would increase the source brightness for 
a given color. We compare the number of h & \ Per 
sources less than 15' away from either cluster core and 
the number between 15' and 25' away. Sources > 15' 
from one cluster and < 15' from the other cluster were 
not counted. 

From the V/V-J diagram, <~ 33% of sources within 
15' of the cluster centers have a CMD position consis- 
tent with cluster membership (~ 767 with J < 14.5). 
In the sample of stars at 15'-25' from the cluster cen- 
ters, <~ 23% lie in a similar zone on the CMD. About 
~ 47% of sources within 15' have a CMD position con- 
sistent with a 13 Myr isochrone (~ 2600 sources with J 
< 15.5) on the J/J-H diagram. Stars near the 13 Myr 
isochrone between 15' and 25' away from either center 
make up a comparable fraction ~ 41%. If we group the 
two populations together <~ 44% of the ~ 11,000 sources 
are of the same age, distance and reddening as h & \ 
Persei. Therefore, within 25' of the cluster centers the 
population of sources with J < 15.5 lying near a 13 Myr 
isochrone is <~ 4700: <~ 2500 closer to h Persei and ~ 
2200 closer to \ Persei. Estimates for the fraction of 13 
Myr sources from the J/J-H diagram are comparable (~ 
27% and 20% for < 15' and 15-25' respectively) to those 
from the V/V-J diagram in the appropriate limit (J < 
14.5): the lower fractional memberships for J < 14.5 are 
likely caused by a high number of foreground K and M 
stars and background M supergiants. Thus, to within 
~ 10-15%, both the V/V-J (J < 14.5) and J/J-H (J < 
14.5, 15.5) CMDs predict comparable estimates of the 
percentage 13 Myr sources at the distance of h & \ Per- 
sci. Furthermore, the percentage of sources lying on the 
13 Myr isochrone in the low-density regions (15'-25') and 
high-density regions (< 15') are consistent within ~ 10% 
of one another. If the 'background' were dominated by 
foreground and background sources then the percentage 
of 13 Myr sources in the high-density regions should be 
much larger than in the low-density regions. 

It is difficult for a random distribution of fore- 
ground/background stars with a wide range of ages and 
spectral types to mimic a ~ 10 yr isochrone, so we con- 
clude that at least some of the 'background' population 
includes young stars with nearly identical ages and dis- 
tances as those in the two clusters. From the poorer 
definition of the cluster along the isochrones in the right 
panels of Figure 8, these stars may be distributed over 
a few hundred parsecs along the line of sight, centered 
roughly on the clusters and probably associated with 
them. The population of sources roughly tracking the 
13Myr isochrone in J/J-H only slowly disappears by ~ 
60' from the cluster centers (not shown), suggesting a 
diameter of about 100 pc. However, while it may be 



possible that the halo population extends to ~ 60' away 
from h & x Per, carefully quantifying its disappearance, 
particularly for the upper main sequence at <~ 13 Myr, 
requires wider angle imaging and deep spectroscopy and 
is beyond the scope of this paper. 

The existence of a halo population of stars at about 
the same age, distance, and reddening as those within 
the cores of h & \ Persei was considered by Schild (1967) 
based on spectroscopy and photometry of the brightest 
sources; S02 also noted that, to V ~ 16, the optical col- 
ors of bright sources beyond 5' of cither cluster center 
and those within 5' of cither cluster center appeared 
quite similar. Our result, probing stars slightly fainter 
than those studied by S02 and over larger spatial area, 
is broadly consistent with both of these references. For 
the rest of the paper, we shall restrict our analysis to 
sources within 25' of the cluster centers. We use the 
V/V-J and J/J-H diagrams to identify sources within 
this radius that appear to have the age, distance, and 
reddening of h & \ Per. 

3.2. JHK S color-color diagrams from 2MASS and 
Mimir 

Now we investigate the near-IR 2MASS colors of all 
sources within 25' of the cluster centers (<~ 11,000) with J 
< 15.5. The J — H / H — K s color-color diagram is shown 
in Figure 9 with the main sequence locus, giant locus, 
classical T Tauri locus, and reddening bands. Stars with 
J-H w 0.0-0.25 and R-K s w 0.0-0.3 will include Be stars 
with near-IR excess from optically-thin gas (Dougherty 
et al. 1991,1994). Excesses around later-type stars with 
J-H > 0.3 and H-K s > 0.35 are more likely to have near- 
IR excess emission from warm dust. The vast majority 
of sources have photospheric colors. Sources off the lo- 
cus are distributed evenly on both sides, suggesting that 
photometric errors may be responsible for any H-K s 'ex- 
cess'. To characterize the population more accurately 
and (later) search for near-IR excess sources at H — K s 
we rely on Mimir data with deeper completeness limits 
at JHK S . 

We analyze Mimir sources with a < 0.2 and J < 16. 
Figure 10 shows the J — H/H — K s color-color diagram 
for all sources in the Mimir survey. The vast majority 
of sources fall along the main sequence locus. Figure 
11 shows data separately for h Persei and \ Persei; h 
Persei has more sources with a slight H — K s excess for 
small J-H colors (< 0.3). The horizontal spread in color 
is noticeably smaller for sources with a < 0.1. Many 
sources have colors consistent with disk excess emission 
(H-K s ~ 0.4 or redder). However, restricting ourselves to 
sources that may be cluster members based on the J/J- 
H diagram, the IR excess population from H-K s through 
J=15.5-16 is extremely small: less than ~ 1%. 

4. IRAC ANALYSIS: 2MASS -[IRAC] COLORS & IRAC-ONLY 

COLORS 

To learn whether excesses are more common at longer 
wavelengths we consider the IRAC data. Restricting our 
analysis to sources within 25' of the cluster centers, we 
have ~ 7000 sources with 5a detections at [3.6] and [4.5]; 
~ 5000 of these have 5a detections at [5.8] and [8.0]. To 
minimize contamination due to large errors, we restrict 
this sample to the sources with if s , [IRAC] < 14.5. 

The R-K s /K s -[3.6] and R-K s / K s -[4.5] diagrams arc 
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shown in Figure 12. In each plot, nearly all sources have 
photospheric colors, with K s -[3.6], [4.5] w 0.0-0.3 or 0.4. 
Sources with if s - [3.6], [4.5] much larger than 0.3-0.4 are 
much redder than normal stellar photospheres. The num- 
ber of 'red' sources is larger at [4.5] than at [3.6] and 
suggests the existence of an IR excess population in h & 
X Per. 

Figure 13 shows the [3.6]-[4.5]/[4.5]-[5.8] and K s - 
[3.6] /-Kg- [5.8] color-color diagrams. In previous studies of 
IRAC colors of pre-main sequence stars (e.g. Hartmann 
et al. 2005), a division between class II and class III T 
Tauri stars occurs roughly at [4.5] -[5.8] <~ 0.2-0.25. The 
h & x Persei color distribution has a red limit of about 
0.3 resulting from the intrinsic dispersion of colors. Many 
sources have colors redder than this limit (if s ,[4.5]-[5.8] 
> 0.3). There are few sources with [4.5]-[5.8] < -0.3. 
The common-baseline i4r s -[3.6]/if s -[5.8] diagram shows 
the very red population even more clearly. Photospheric 
sources have if s -[5.8] < 0.3-0.4. The plot also suggests 
that there may be some sources with red colors at K s - 
[5.8] but not if s -[4.5], judging from the relatively larger 
number of sources. 

The distribution at [5. 8] -[8] (Figure 14, left panel) con- 
tinues to show a very gradual, as opposed to abrupt, 
transition in colors from photospheric ([5.8]-[8] < 0.2) to 
very red ([5.8]- [8] > 0.4, see Allen et al. 2004). In Hart- 
mann et al. (2005) there are very few sources with [5.8]- 
[8]=0. 25-0.4. The lack of any gap in h & X Per is likely 
due to larger photometric errors at [8]. The same plot 
for 10er sources (Figure 14, right panel) shows a much 
smaller dispersion in [5.8]- [8] colors containing sources 
with [5.8]-[8]=0.25-0.4. There is also a substantial num- 
ber of 10cr detections with [5.8]-[8] > 0.4. 

Constructing colors from the 2MASS K s band and 
IRAC bands at [5.8] and [8] provides the clearest evi- 
dence for an IR excess population in h & x Persei. Figure 
15 shows the K s -[3.6]/ K s -[8] color-color diagram, which 
is our closest analogue to the K-L/K-N diagram used by 
Kenyon & Hartmann (1995) to distinguish Class II and 
III T Tauri stars in Taurus. Photospheric sources appear 
to have K S -[S}< 0.4; with K s -[8] > 0.4 there is a clear 
population of IR excess sources. The ^-[5.8]/^- [8] di- 
agram also shows an excess population for if s -[5.8,8] > 
0.4. 

The population of very red, IR excess sources in K s - 
[IRAC] is statistically significant. The mean and disper- 
sion in K s -[4.5] arc - -0.02 ± 0.06; 0.07 ± 0.1; and 0.07 
± 0.1 at J=12, 14.5, and 15. For if s -[8] color these values 
are 0.06 ± 0.1, 0.15 ± 0.18, and 0.24 ± 0.19; the values 
for if s -[3.6] are similar to AT s -[4.5] and those for if s -[5.8] 
are in between if s -[4.5] and [8]. The median colors are 
nearly identical in all cases. In all if s -[IRAC] colors there 
exists a substantial population of sources that are more 
than 2-5<7 redder than the mean color while a correspond- 
ing blue population does not exist. It is clear from the 
IRAC colors that h & x Persei harbors a significant IR 
excess population. 

5. ANALYSIS OF THE IR EXCESS POPULATION 

Now we quantify the fraction of sources with excesses 
from 2MASS and IRAC photometry and thus the pop- 
ulation of sources with circumstellar disks. With no ev- 
idence for IR excess at K s , we estimate the fraction of 
sources with IR excess at [4.5], [5.8], and [8], using K s as 



a common short wavelength baseline and then relate the 
IR excess population to intrinsic stellar properties via the 
2MASS J band. The J filter should have emission dom- 
inated by the stellar photosphere, especially for sources 
older than 10 Myr that typically are not actively accret- 
ing. Converting from J magnitudes to stellar properties 
is also more straightforward than with optical filters be- 
cause J is less affected by reddening. We analyze the 
IR excess population as a function of J in 0.5 magnitude 
bins. 

In Section 5.1, we describe our two methods for select- 
ing sources as h & x Persei members. Next, we describe 
our criteria for identifying a source as an IR excess source 
in section 5.2. We estimate the size of the IR excess pop- 
ulation using a highly restrictive model for membership 
and IR excess identification and then using a less restric- 
tive model. We describe both models in section 5.3. 

The estimates from both the restrictive and less restric- 
tive models show excellent agreement. The relative size 
of the IR excess population to the total population in- 
creases with J magnitude. Most of the IR excess sources 
have J=14-15; very few excess sources are brighter than 
J=13.5-14. This result implies that the frequency of disks 
around stars with ages 10-15 Myr is likely higher for lower 
stellar masses. The IR excess population is consistently 
larger at longer IRAC wavelengths than at shorter wave- 
lengths: this behavior is expected if circumstellar disks 
clear from the inside out. At least ~ 4 — 8% of sources 
with J=14-15 (~ 2.2 — 1.4M Q ) have IR excess indicative 
of a circumstellar disk. We also identify sources that 
have IR excess at both [5.8] and [8]. About - 2 - 3% of 
sources have strong IR excess at both long wavelength 
IRAC channels for sources as faint as J=15.0. 

5.1. Sample Selection 

Our first task is to remove probable non members 
of h & x Persei. After removing sources identified as 
AGN/aromatic-emission galaxies we used the V/V-J and 
J/J-H color-magnitudes to remove sources inconsistent 
with the 13 Myr isochrone (Section 2.3). The most re- 
strictive, cautious approach is to use only the sources 
with optical photometry from S02 because the long- 
baseline V-J color is better at identifying non members 
than the shorter-baseline J-H color. For this approach 
we require sources lying within 0.3 magnitudes (~ 5a 
errors added in quadrature) of the isochrone in V 

\V{source) — V(isochrone)\ < 0.3. (1) 

Thus our sources consistent with h & x Per membership 
should fall within a band 0.3 magnitudes brighter and 
fainter than the nominal reddened isochrone. However, 
the V magnitude data from S02 is complete only to ~ 
14.5 in J. Thus, restricting our analysis to optical sources 
eliminates stars with J > 14.5 or M < 1.6 M©. 

A less restrictive approach uses the 2MASS J-H color 
to identify ~ 13 Myr, 2.34 kpc distant stars. Sources 
in unresolved binary systems should be < 0.75 magni- 
tudes too bright for their color. The faint limit, however, 
should be unchanged. Our criteria are then: 

J(source) — J {isochrone) > —0.75, < 0.3. (2) 

This 2MASS sample is complete to J=15.5, which in- 
cludes stars with M > 1.3M Q . 
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5.2. IR Excess Criteria 

A color threshold is set for a source to be considered 
an IR excess candidate using the 2MASS - [IRAC] color- 
color diagrams and Kenyon & Hartmann (1995) photom- 
etry and color table guide. For the [4.5] channel 
we use the K, 
nel we use K, 



3.6]/K 
3.6]/AV 



4.5] diagram, for the [5.8] chan- 
5.8], and for the [8] channel we 
use K s -[3.6]/K s -[8] (Figures 11, 12, and 14). The main 
distribution of sources in these diagrams has K S -[IRAC] 
< 0.4. For comparison, Kenyon & Hartmann found that 
sources with K-L < 0.3 — 0.4 were typically class III T 
Tauri stars with weak H a emission and little evidence 
of circumstellar dust at L band (see also McCabe et al. 
2006, figures 1 and 2). According to the Siess et al. 
(2000) isochrones, 13 Myr sources more massive than ~ 
IMq should also not have K-L> 0.2 — 0.25, assuming the 
extinction derived by BK05. Sources with K s - [IRAC] 
> 0.4 are then likely IR excess source candidates. 

To avoid defining sources with large photometric errors 
as excess sources, we make the additional requirement 
that a source's color must be redder than the absolute 
threshold plus the source's photometric error, a, in the 
IRAC bands, or 

K s - [IRAC] > 0.4 + a(IRAC). (3) 

Fainter IRAC sources (~ 14-14.5) then typically must 
have .ff s - [IRAC] > 0.5-0.6 to be classified as excess 
sources. As a sanity check, we also require that the next 
longest wavelength IRAC band show at least a marginal 
excess, K S -[IRAC] > 0.3, if a 5a detection was made at 
such a channel. 

Finally, we remove potential completeness-related bias 
and calculate errors on the size estimates of the IR ex- 
cess population. First, because the completeness limits 
from J through [8] vary, we apply a uniform cutoff across 
2MASS and IRAC bands and restrict ourselves to sources 
with if s , [IRAC] < 14.5, as the count rate for 5a detec- 
tions at [8] falls to half-peak values by 14.5. Because 
sources with K s ~ 14.5 have 3-K s ~ 0.4-0.7, we add a 
further restriction that sources must have J< 15.0. Es- 
timates for the size of the IR excess population in each 
IRAC band are then drawn from a single population. 
Nominally, the errors for the IR excess fraction in each 
0.5 J magnitude bin are calculated from Poisson statis- 
tics, <7<~ \/#/N, where # is the number of sources with 
IR excess and N is the number of sources in the bin with 
detections in 2MASS and at a given IRAC band. We re- 
tain this error estimate as an upper bound on the IR ex- 
cess fraction. However, there may be sources in our cov- 
erage that are detected in 2MASS but not in IRAC and 
are brighter than our J magnitude cutoff. Thus, these 
sources would not be included in our calculation for N 
because of their non-detection in IRAC. This could bias 
our results in favor of detecting more IR excess sources 
and fewer photospheric sources near magnitude 14.5 and 
thus is another source of error. 

To put a limit on this error, we calculate the num- 
ber of sources in each bin detected at both J and K s . 
If the IRAC sample is 100% complete, this number M 
should be equal to N, otherwise it will be larger. We 
then assume that any source in M not detected in N is 
not an IR excess source and divide the total number of 
IR excess sources by M. This procedure yields an abso- 
lute lower limit on the fraction of sources with IR excess 



in each 0.5 magnitude bin. We subtract this value from 
the nominal fraction to produce an error estimate due to 
sample completeness. The larger of the two errors, the 
Poisson error and the completeness error, is chosen as 
the lower bound on the IR excess fraction. 

5.3. Two Models for Analyzing the IR Excess 

Population 

Equipped with our two approaches for sample selec- 
tion, our method for identifying an IR excess source, 
and our brightness cutoffs in 2MASS & IRAC, we now 
describe our two models for analyzing the IR excess 
population. The models are summarized in Table 3. 

Model 1 (Restrictive) - We require that sources 
have optical counterparts from S02 and that they are 
within 0.3 magnitudes of the 13 Myr V/V-J isochrone 
(Equation 1). We require that IR excess candidates 
fulfill the condition set in Equation 3. We investigate 
the IR excess population at [4.5], [5.8], and [8] for 
sources with J < 14.5. Errors for the size of the IR 
excess population in all models assume Poisson statistics. 

Model 2 (Less Restrictive) - We require that sources 
have near-IR colors consistent with the 13 Myr J/J-H 
isochrone as described by Equation 2. We require 
that the IR excess candidates fulfill the condition set 
in Equation 3. The IR excess population is analyzed 
through J=15. 

As additional checks of the reliability of the data, we 
examined the individual frames for the sources with ap- 
parent excesses and considered as highest weight the ex- 
cesses detected at both [5.8] and [8]. We quantify the 
excess population at both [5.8] and [8] as well. 

5.4. IR Excess Population: Restrictive Model 

We show the size of the IR excess population relative to 
the total population as a function of J magnitude for our 
restricted model (Model 1) for K s -[4.5], K s -[5.8], and K s - 
[8] in Figure 16. The most striking feature in the diagram 
is the complete lack of IR excess sources at any band for 
J < 13.5: 0/250 sources with J=ll-13.5 have IR excess. 
We only begin to detect excess sources at [8] by J=13.9. 
The IR excess population at [4.5] and [5.8] only starts to 
appear by J ~ 14.05. While our approach may under- 
estimate the fraction of IR excess sources near our faint 
limit of J=14.5 (M ~ 1.6M ), sources with J=ll-13.5 
typically have small photometric errors (< 0.05 mag) at 
all 2MASS and IRAC bands. Thus, we have not missed 
IR excess sources brighter than J=13.5 (M ~ 2.7M Q ) 
from our sample. 

We note two other features in the IR excess popula- 
tion. The [8] IR excess fraction is consistently larger 
than at [4.5] and [5.8]. This result is consistent with 
an inside-out clearing of circumstellar disks. Stars with 
J=14-14.5 (2.2-1.6 M Q also tend to have a larger frac- 
tion of IR excess sources than stars with J=13.5-14. For 
J=14-14.5, 14/409 (3.4 +0.9,-1%) sources have IR excess 
at [8]; for J=13.5-14 only 2/169 (1.2 ± 0.8%) of sources 
have [8] excess. This result coupled with the lack of ex- 
cess sources at J< 13.5 suggests that the frequency of IR 
excess sources is related to stellar properties. 

5.5. IR Excess Population: Less Restrictive Model 



h & x Persei 



9 



Now we investigate the IR excess population as a func- 
tion of J magnitude for our less restrictive model (Model 
2; Figure 17). While this approach may introduce more 
non members, the larger number of cluster members in 
this sample (to J< 15.5) yields a larger statistical sig- 
nificance for our results. There are no IR excess sources 
for J< 13.5: 0/399 with J=ll-13.5 have IR excess at [8] 
or at shorter wavelengths. Thus the lack of IR excess 
sources at IRAC bands to J=13.5 is model independent. 

The IR excess population from this sample appears 
remarkably similar to the restrictive sample where the 
two overlap (to J=14.5). In both models the [8] excess 
population is consistently larger than that at [4.5] and 
[5.8]. The IR excess population at [8] increases from 
1.6 ± 0.6% (8/496) of the total population from J=13.5- 

14 to 3.5 +0.6,-0.9 % (34/983) at J=14-14.5. This result 
agrees with the estimates (complete to J=14.5) derived 
in section 5.4. Thus to J— 14.5, our estimates of the IR 
excess population's size is not sensitive to whether we 
use the long-baseline V/V-J or the short-baseline J/J-H 
diagram to constrain cluster membership. For J=14.5- 

15 the IR excess population at [8] is <~ 8.1 +1.1,-4.5% 
(50/618). The population of sources with J=14.5-15, 
K s < 14.5, and [8] detected with errors < 0.2 is smaller 
than the total population of sources with J=14.5-15 and 
K s < 14.5: this accounts for the relatively large lower 
bound error. We note that if the sample were restricted 
to -K" s ,[8] < 14.25, where completeness is better, the per- 
centage of excess sources is roughly halfway between our 
measured value of 8.1% and the lower limit of ~ 4% at 
~ 6.3%. Thus, the IR excess population at [8] near the 
brightness limit of our analysis comprises about 4-8% of 
the total population. 

The excess population at X s -[4.5] is consistently the 
smallest, ranging from zero through J=13.5 to 0.2 
±0.2%, 1.4+0.4% (16/1117), to 1.0 +0.5,-0.6% (11/1138) 
of the total population from J=13.5-14, 14-14.5, and 
14.5-15 respectively. Over the same range in J magni- 
tude the [5.8] excess population represents respectively 
zero, 0.6+0.3%, 1.6+0.4% (19/1214), and 3.3 +0.5,-0.6% 
(38/1150) of the total population. 

5.6. Sources with Clear IR Excess at both [5.8] and [8] 

If the IR excess population is dominated by sources 
with optically thick circumstellar disks extending to the 
magnetospheric truncation radius (Kenyon et al. 1996), 
then the IR excess population at [5.8] and [8] should 
be nearly identical. However, in both models for ana- 
lyzing the IR excess population, the fraction of sources 
with IR excess at [8] is consistently larger. It is plausible 
that sources may not have IR excess at short-wavelength 
IRAC bands but have excess beyond <~ 5 — 10/zm, espe- 
cially at 10-20 Myr. For example, TW Hya, a 10 Myr T 
Tauri star, has photosphcric colors through [4.5], and a 
clear IR excess only at [8] (Hartmann et al. 2005). Other 
evolved T Tauri stars, such as the 'transitional' T Tauri 
stars in Taurus (Kenyon & Hartmann 1995), have SEDs 
consist with a star+circumstellar disk where the disk has 
an inner hole. Sources with strong IR excess longward of 
[5.8] may also be explained by debris produced by planet 
formation (Kenyon & Bromley 2004; Currie et al. 2006). 
However, unlike for TW Hya and similar systems, we do 
not have constraints beyond 8[im on IR excess candi- 
dates in h & x Per since the MIPS completeness limit is 



typically too bright to detect these candidates (Currie et 
al. 2006; Balog & Currie et al., in prep.). 

Therefore sources with both [5.8] and [8] excess set an 
absolute lower limit on the size of the IR excess popula- 
tion. We set the excess criteria to K s -[5.8,8] > 0.4 and 
use our sample from Model 2. The IR excess popula- 
tion at both [5.8] and [8] begins at J=13.8 and represents 
<~ 1% of the population for sources with J=14-14.5. This 
percentage increases to <~ 2.2% for J=14.5-15, slightly 
less than but still consistent with the excess fraction at 
[5.8] alone from Model 2. We also recover the same IR 
excess population vs. J magnitude relation in models 1 
and 2, albeit at a lower statistical significance. 

Thus, for sources with M <~ 1.4 — 1.6M Q at least ~ 
2 — 3% have IR excess at both [5.8] and [8]. If we include 
sources with photospheric colors at all wavelengths short 
of 8 fxm then about <~ 4 — 8% of sources with J=14.5-15 
have IR excess. 

5.7. Distribution of IR excess sources on V/V-J and 
J/J-H Color-Magnitude Diagrams 

Here we consider estimates of the IR excess population 
without cluster membership criteria. If IR excess sources 
are distributed randomly in V/V-J or J/J-H space our 
IR excess population estimates might be invalid. 

Figure 18 shows our results. In the top panel of Figure 
18 IR excess sources with V band photometry from S02 
preferentially lie along the 13 Myr isochrone. They are 
far more clustered around the isochrone than the distri- 
bution of all sources. This trend also prevails for sources 
without V band photometry as shown in the bottom 
panel of Figure 18. Here, the IR excess population is pref- 
erentially clustered along the 13 Myr J/J-H isochrone. 
Just over half of the sources within 25' of h & % Persei are 
foreground/background stars, so the sizes of the h & x 
Per population and foreground/background populations 
are comparable. However, the IR excess population for h 
& X Per cluster and halo population members is clearly 
larger than the population of non members. This result 
is expected if our membership and IR excess criteria arc 
valid. 

5.8. Statistical Trends in the IR Excess Population of h 
& x Persei 

To conclude, analyzing the excess populations at K s 
-[4.5], [5.8], and [8] according to our two models leads to 
the following trends. First, both models show a complete 
lack of h & x Persei sources with J< 13.5 and IR excess 
at any band. This result suggests that almost all inner 
disks (< about 1 AU from the star) around massive stars 
(> 2.5 — 3Mq) disappear prior to ~ 10 Myr. Second, 
both models show a general increase in the IR excess 
population with J magnitude. This suggests that a small 
number of disks around lower mass stars (to ~ 1.4 — 
I.QMq) have lifetimes longer than <~ 10 Myr. Third, to 
J=14.5 (<~ 1.6M Q ) about 3-4 % of sources have IR excess 
in the IRAC bands. 

The deeper completeness and stronger statistical con- 
straints from Model 2 separately yield the following ad- 
ditional results. The IR excess population is larger at [8] 
than at [5.8] and, especially, [4.5]. This result is consis- 
tent with an inside-out clearing of disks. About 4-8% of 
sources slightly more massive than the Sun (J=14.5-15, 
<~ 1.4 — 1.6M Q ) have IR excess. Tighter constraints on 
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the cluster membership from optical spectroscopy and 
deeper IRAC and MIPS observations will be required to 
more precisely determine the IR excess population. 

6. SUMMARY, FUTURE WORK, & DISCUSSION 

We have conducted the first deep IR survey of the dou- 
ble cluster, h & \ P er - The survey combines 2MASS 
JHK S data with deeper JHif s data from the Lowell 
Mimir camera and Spitzer IRAC 3.6-8fim data cover- 
ing a region with an area of roughly 0.75 deg 2 . To derive 
the fraction of stars with IR excesses, wc considered the 
observed stellar surface density distribution, optical and 
infrared color-magnitude diagrams, and infrared color- 
color diagrams. This analysis yields the following robust 
results: 

• For stars with J < 15.5, the stellar surface density 
peaks at 15 arcmirC 2 . In both clusters, the den- 
sity falls to background levels of ~ 2.7 arcmin~ 2 
<~ 15' away from the center peak. About ~ 47% 
of the stars within 15' are cluster members. Com- 
parisons of color-magnitude diagrams show that ~ 
41% of the stars at 15-25' from the cluster centers 
are of approximately the same age, distance, and 
reddening. In both cases, the Siess et al (2000) 
evolutionary tracks suggest ages of <~ 13 Myr for 
cluster stars. 

• For stars with J < 13.5, the two clusters have es- 
sentially no stars with IR excess at A < 8 /im. This 
result demonstrates that disk emission with R < 1 
AU disappears from stars with M > 2.5-3 M Q on 
timescales of < 10-15 Myr. 

• Many stars with J =13.5-15 have IR excesses: > 
2-3% of the stars in this magnitude range have 
excesses at both [5.8] and [8.0]. Combined with 
the lack of detectable excesses for brighter stars, 
these results suggest that the disk evolution time 
depends on the stellar mass, with lower mass stars 
having longer disk evolution timescales. 

In addition to these firm conclusions, our analysis indi- 
cates the following more tentative results: 

• If the requirement that sources have excesses at 
both [5.8] and [8] is relaxed, the fraction of sources 
with IR excesses increases from 3-4% at J = 14-14.5 
to 4-8% at J=14.5-15. 

• At a fixed brightness, a larger fraction of sources 
appear to have excesses at longer wavelengths, 5.8- 
8 /zm, than at shorter wavelengths, 4.5 /im. This 
result is consistent with an inside-out clearing of 
protoplanetary disks as proposed in standard the- 
ories of planet formation. 

These results demonstrate that h & \ P er are an ex ~ 
ccllcnt laboratory for studying the evolution of circum- 
stellar disks at 10-15 Myr. With > 5000 sources more 



massive than the Sun and a potentially significant halo 
population outside the main cluster boundaries, h & \ 
Per are at least as populous as the Orion Nebula Clus- 
ter (ONC; Hillenbrand 1997) and are the most populous 
nearby clusters in this age range within 2-3 kpc. De- 
tailed comparisons of these clusters with the ONC and 
other young clusters provide robust tests of the univer- 
sality of the IMF (S02) and the evolution of small scale 
structures within the clusters (Hillenbrand & Hartmann 
1998, BK05). 

Further progress on understanding the nature of the 
IR excess population requires deeper IR imaging surveys 
and comprehensive spectroscopic surveys. Current near- 
IR imagers on 4-10 m class telescopes can reach 0.5 Mq 
stars (J ~ 17.5) in a few nights. When combined with a 
deeper IRAC survey ([3.6, 4.5] - 16 - 17, [5.8, 8] - 16), 
these data would yield measures of the IR excess for low 
mass stars where the predicted timescales and outcomes 
for disk evolution and planet formation are much dif- 
ferent than at 1-3 M Q (Plavchan et al. 2005; Laughlin 
et al. 2004). Modern multi-object optical spectrgraphs 
routinely acquire high S/N spectra of V = 19-20 stars in 
30-45 min. With current spectroscopic samples complete 
only to V ~ 15-16 (S02, BK05), a deep spectroscopic sur- 
vey enables more reliable measures of cluster membership 
and deep searches for Ha emission for stars with masses 
of - 0.5-2 M fl . 
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TABLE 1 

MlMIR OBSERVATIONS COVERAGE 



Field 


Co 


So 


Date (2005) 


I 




57°8'35" 


INUvLlllULI '-tO 


2 


9h1 Qrn1 6 8*i 


O 1 o o oo 




3 


9h1 Qrn^7 3<3 


57°8'35" 


IN (J Vt^IllDUl '-t U 


4 




57°8'35" 


IN UVCLIIDL-I U 


5 


9h1 Qrn1 6 8<s 


57°5'35" 




6 


2hl9m37.2s 


57°5'35" 


November 6 


7 


2hl9m57.6s 


57°5'35" 


November 6 


15 


2h22m4.3s 


57°8'35" 


November 30 


16 


2h21m43.2s 


57°8'35" 


December 5 


17 


2h21m22.1s 


57°8'35" 


December 1 


19 


2h20m39.7s 


57°8'35" 


December 5 


20 


2h20ml8.5s 


57°8'35" 


December 5 


22 


2h22m4.3s 


57°5'35" 


December 1 


23 


2h21m43.9s 


57°5'35" 


December 1 


24 


2h21m23.5s 


57°5'35" 


December 4 


25 


2hl9ml7.0s 


57°10'0" 


December 5 



TABLE 2 
H & x Persei data 



a 


5 


J 


H 


K s 


[3.6] 


[4.5] 


[5.8] 


[8] 


<r(J) 


<r(H) 




<r([3.6]) 


<r([4.5]) 


<r([5.8]) 


*([»]) 


34.9281 


56.9391 


14.394 


13.875 


13.691 


13.613 


13.619 


13.451 


13.583 


.034 


.038 


.049 


.003 


.428 


.07 


.118 


35.0360 


56.8215 


14.229 


13.654 


13.656 


13.517 


13.562 


13.883 


99.0 


.035 


.042 


.048 


.039 


.056 


.098 


99.0 


35.0483 


57.0313 


14.516 


13.983 


13.857 


13.752 


13.781 


13.771 


13.563 


.034 


.037 


.057 


.006 


.017 


.104 


.124 


35.0523 


56.7621 


12.421 


12.001 


11.922 


11.816 


11.844 


11.920 


11.903 


.024 


.030 


.023 


.057 


.005 


.022 


.030 


35.1106 


56.8819 


8.795 


8.049 


7.836 


8.324 


7.9730 


7.787 


7.731 


.018 


.024 


.018 


.014 


.168 


.002 


.003 



NOTE. — First five entries in our photometry catalogue from 2MASS and IRAC 



TABLE 3 

Models for estimating the IR excess population 

Model Membership Determination Membership Cutoff J cutoff 

1 (more restrictive) V/V-J colors 1 |V-V(isochronc)| < 0.3 14.5 

2 (less restrictive) J/J-H colors 2 J-J(isochrone) > -0.75, < 0.3 15.0 



NOTE. — Brief description of the two models for identifying IR. excess sources. In the first model we require that sources have optical data from S02 and use the V/V-J 
diagram to constrain cluster membership. We analyze sources to .1=14.5. In the second model we use the J/J-H diagram from 2MASS to constrain cluster membership. We 
analyze sources to J=15.. The membership cutoff refers to how far away a source can be from the 13 Myr isochrone and still be classified as an h & x Persei member. The 
cutoff limits in the second model include the effect of binarity. 

a Data from S02 and 2MASS 
^Data from 2MASS 
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Fig. 1. — Coverage map for Mimir (cyan boxes) and IRAC (purple dash-dotted line) observations plotted over 2MASS sources in h &c x 
Per. The center of h Persei (2 h 22 m 4.3 s 57°8'25") is in the extreme right cyan box and the x Persei center (2 /l 22 m 4.3 s , 57°8'35") is in the 
two leftmost boxes. The areas covered by Mimir and IRAC are ~ 144 arc-min 2 and ~ 0.75 deg 2 respectively. 
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Fig. 2. — J (top) and K s (bottom) magnitude and error distributions for 2MASS sources within 25' of the h & \ Per centers. The dotted 
lines are for sources closer to h Persei and the dash-dotted lines are sources closer to \ Persei. The vertical dashed lines are the published 
lOcr sensitivity limits for 2MASS: J=15.8, K S =IA.3. The small turnover near J ~ 14.5, K a ~ 14.25 results from the large number of sources 
near the centers of h & \ Per that follow the ~ 13 Myr isochrone. The total sample is complete to J=15.5 and K s ^ 15.0 (~ 1.3Mq). falls 
to half the peak value by J ~ 15.75 — 16 and K s ^ 15.25. Errors at J=15.5 and K a =15.0 are less than ~ 0.1 and ~ 0.12 respectively. The 
distribution in K s band shows reasonable agreement with the published 10<r limits while those in J fall slightly short but are reasonable 
through J-~ 15.5. We restrict the data analysis to sources with J < 15.5. 





15 ^ 16 17 15 14 15 16 



Fig. 3. — (top-left) J distribution for h Persei sources (dotted line), \ Persei sources (dot-dashed), and all h & x P er sources (solid 
line) from Mimir . There is a small enhancement of sources near J ~ 14.0 — 14.5 which are cluster sources on a 13 Myr isochrone. The 
number counts in J peak at ~ 16.0 — 16.25 which corresponds to stars with 1.1 — 1.2Mq at the adopted distance of h & \ Per (~ 2.34 
kpc). However, the distribution does not fall to half its peak value until J ~ 16.75. We see a secondary peak at J ~ 14.5. (top-right) 
K s distribution for h & x Persei sources. This distribution has peaks at 14.25 and 15.5. (bottom) The 2MASS J and K s magnitudes vs. 
Mimir J and K s magnitudes for faint 2MASS sources. Dotted lines show ~ 0.2 magnitude deviations. The vast majority of sources fall 
within the dotted lines; the majority of sources falling outside the lines have fainter Mimir magnitudes. This is consistent with resolving 
binary stars, a likelihood given Mimir's higher spatial resolution. We restrict our analysis of the Mimir data to sources with J < 16. 
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Fig. 4. — IRAC Mosaic image of the h & \ Persei region at 3.6 fim. The centers of h Persei and \ Persei are in the center right and 
center-left portions of the image respectively. The contrast is set to 98%. Diffraction spikes are visible on several of the brightest sources. 
The Spitzcr/IRAC coverage is ~ 1 square degree on the sky. 
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Fig. 5. — Distributions of 5<r detections at [3.6], [4.5], [5.8], and [8]; error distributions at [4.5] and [8]. The [3.6] and [4.5] data are 
complete to ~ 14.5 but have a substantial population to m(3.6,4.5)~ 15. The [5.8] data are complete to ~ 14.25 with a substantial 
population to ~ 14.5. The [8] data are complete to ~ 13.75 and falls below half the peak value by ~ 14.5. The vast majority of sources in 
[4.5] have errors < 0.15 to m(4.5)=15; sources in [8] have errors < 0.2 to m(8)=14.5. 
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Fig. 6. — The completeness in each IRAC band as a function of 2MASS J magnitude (to J=15.0). Due to the small matching radius 
used to combine the 2MASS and IRAC data sets, the completeness level for IRAC oscillates between 90 and 100% for J< 14. Except for 
the [8] band for sources with J=14.5-15, all IRAC bands are better than 85% complete. In the IRAC analysis section (Sect. 5) wc account 
for completeness errors in estimating the IR excess population. 
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Fig. 7. — (Top) Spatial plot of the star surface density for sources with J < 15.5 from to 90% in increments of 10%. The centers of 
both h & x Persei (right and left peaks, respectively) arc clearly visible and are separated by about ~ 26'. The clusters are ~ 20'-30' 
across with substantial asymmetries, most notably close to the \ Persei core at slightly larger declinations. (Bottom) Radial density plot 
of sources with J < 15.5. The counts reach the median background density (~ 2.7 arcmin~ 2 ; dash-three dots) by ~ 15 — 25' away from 
both h Per and \ Per. 
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Fig. 8. — V/V-J (top) and J/J-H (bottom) color-magnitude diagrams of the 2MASS and S02 sources with <r < 0.2 for sources within 
15' (left) and those between 15' and 25' (right) from the h & \ Per cluster centers. Overplotted are the Siess et al. (2001) isochrones 
for a 13 Myr cluster (dash-three dots) with a 'reddened' isochrone (E(B-V)~ 0.52; solid line) and appropriate reddening vectors (Ay ~ 
1.4 E(V-J); Aj ~ 2.49 E(J-H)). For the J/J-H diagram, the reddened pre-main sequence tracks for 0.8-7 Mq stars from Bernasconi et al. 
(1996, dot-dashed) and those for lower mass stars from Baraffe et al. (1998, long dashes) are also plotted. Sources between 15' and 25' 
away from the cluster centers appear to contain a population tracking the isochrone by V> 14.5 and J> 13.5. 
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Fig. 9.— J-H vs. H - K a colors from 2MASS for 11048 sources with J < 15.5 . Overplotted are the classical T Tauri (Meyer et al. 1997) 
and giant loci (Kenyon & Hartmann 1995). Crosses represent the change in colors that a source with photospheric J-H/H-fsfs colors of 
~ 0.55/0.1 would have if a 1500 K blackbody (from a circumstcllar disk) contributed 5%, 25%, and 50% of the total flux at 3.6 fim. 
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-10 12 

H-K s 

Fig. 10. — J-H vs. H - K s colors from Mimir complete for 657 stars with to J=16 (~ 1.1M©, d~2.34 kpc) . Filled diamonds represent 
sources detected at the 10<r (<r < 0.1) level in all bands, triangles represent sources with errors less than a = 0.2 in the three bands. 
Overplotted are the classical T Tauri (Meyer et al. 1997) and giant (Kenyon & Hartmann 1995) loci. 
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Fig. 11. — J-H vs. H - K s colors from Mimir for h Per (left) and \ Per (right). The symbols are the same as in previous plots. 
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Fig. 12.— The H-K s /K s -[3.6] color-color diagram (left), and H-K s /K s -[4.5] diagram (right). The distribution of H-[3.6], H-[4.5], 
K s -[3.6], and K s -[4.5] colors show a potential IR excess population with i^ s -[IRAC] > 0.3-0.4. 
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Fig. 13.— [3.6]-[4.5]/[4.5]-[5.8] (left panel) and K a -[3.6]/K a -[5.8] color-color diagrams (right panel) of sources with J > 11, < 15.5 and 
m(5.8)< 14.5. In the bottom diagram the dotted lines denote sources with [3.6]- [5.8] =0-0.3. The photosphcric population appears to have 
a red edge [4.5]-[5.8]=0.3 and K S -[5.8]=0A. In Section 5 we count sources with A" s -[5.8]> 0.5 as IR excess candidates. 
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Fig. 14. — (Left panel) Four-channel IRAC colors for sources within 25' of h & \ P er centers. The population of sources in the clear 
'break' between Class II and III T Tauri stars from 0.25-0.4 is large compared to the overall population of sources with IR excess (5.8-8 > 
0.4). (Right panel) The same four-channel IRAC color distribution among sources with small (cr < 0.1) errors in [8]. There are sources in 
the 'break' between Class II and Class III colors from [5.8]-[8]=0.25-0.4 and excess sources beyond [5.8]-[8]~ 0.4. 
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Fig. 15.— K s -[3.6]/K 3 -[8] (left panel) and K a -[5.8]/K s -[8] (right) color-color diagrams for sources with J=ll-15.5 and [8]< 14.5. The 
dotted lines in the bottom panel bound sources with [5.8]-[8]= 0-0.4. In both plots there is a red edge to the main population at i<T s -[8]~ 
0.4 and many IR excess sources with K s -[8]> 0.4. 
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Fig. 16. — The disk fraction at 4.5/im, 5.8/im, and 8/um as a function of J magnitude according to Model 1. Fractional excesses are 
estimated in 0.5 magnitude bins from 10.5 through 14.5, which corresponds to a range in mass of 7-1.6 Mq for d~ 2.34 kpc. The IR-excess 
population is larger for fainter sources. We detect no excess sources with J< 13.5. About 3.5% of sources ~ 1.6 Mq have IR excess and 
may harbor circumstcllar disks. 
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Fig. 17. — The disk fraction at 4.5/im, 5.8fim, and 8 fim as a function of J magnitude for Model 2. Fractional excesses are estimated 
in 0.5 magnitude bins corresponding to a mass range of 7-1.4 Mq. The results are consistent with those in Figure 16 (Model 1). The IR 
excess population is clearly larger for fainter sources to J=15.0. There are no IR excess sources for J< 13.5 (~ 2.7Mq), which implies that 
disk lifetimes are larger for lower mass stars (to ~ 1.4Mq). The IR excess population is also consistently larger at progressively longer 
wavelengths, as expected if the disk lifetime in the hotter, inner disk is shorter than in the cooler, outer disk. About 4-8% of sources to 
J=15 have IR excess and may possess circumstcllar disks. 
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Fig. 18. — The distribution of IR excess sources on color-magnitude diagrams of all sources. Both plots show the i\\,-[4.5] (filled 
diamonds), X s -[5.8] (triangles), and i\\,-[8.0] (squares) excess and J<15.5. Arrows correspond to the appropriate reddening vectors as 
in Figure 8. In the top panel, the V/V-J distribution for S02 sources with the IR excess sources is overplotted. The reddened 13 Myr 
isochrone (dash-three dots), the isochrone for unresolved binaries (top solid line) and the unreddened isochrone (bottom solid line) are 
shown, assuming Ay^3. 12 E(B-V)~ 0.16 and E(V-J)~ 1.16 (Bessel & Brett 1988). The vast majority of IR excess sources are consistent 
with the 13 Myr isochrone. The bottom panel shows the 2MASS J/J-H diagram of for IR excess sources with . The reddened 13 Myr 
isochrone (assuming Aj ~ 0.46, E(J — H) ~ 0.185 extinction) the 13 Myr isochrone for unresolved binaries are shown as the bottom and 
top solid black lines, respectively. The unreddened 13 Myr isochrone is also shown (dash-three dots). Again the vast majority of IR excess 
sources are consistent with the 13 Myr isochrone and thus arc consistent with cluster membership. 



